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ABSTRACT 

Wc describe the construction of a suite of galaxy cluster mock catalogues from N-body 
simulations, based on the properties of the new ROSAT-ESO Flux-Limited X-Ray 
(REFLEX II) galaxy cluster catalogue. Our procedure is based on the measurements 
of the cluster abundance, and involves the calibration of the underlying scaling rela- 
tion linking the mass of dark matter haloes to the cluster X-ray luminosity determined 
in the ROSAT energy band 0.1 — 2.4 keV. In order to reproduce the observed abun- 
dance in the luminosity range probed by the REFLEX II X-ray luminosity function 
(0.01 < L x /{10 4i eigs- 1 h- 2 ) < 10), a mass-X ray luminosity relation deviating from 
a simple power law is required. We discuss the dependence of the calibration of this 
scaling relation on the X-ray luminosity and the definition of halo masses and analyse 
the one- and two-point statistical properties of the mock catalogues. Our set of mock 
catalogues provides samples with self-calibrated scaling relations of galaxy clusters to- 
gether with inherent properties of flux-limited surveys. This makes them a useful tool 
to explore different systematic effects and statistical methods involved in constrain- 
ing both astrophysical and cosmological information from present and future galaxy 
cluster surveys. 
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1 INTRODUCTION 

Mock catalogues have become a key tool to assess the sta- 
tistical methods employed in the analysis of redshift galaxy 
surveys, mainly concerni ng the study of the large-scale struc- 
ture of the Un i verse (e. 



haloe s (e.g IWarren et al.ll20o4 iTinker Il200sl : ICrocce et all 
2009) as well as on the non-linear evolution of gravitational 



lg tnc study 01 trie large-scale struc- 
Cole et alJboollCabre et al.ll2009l: 



Evrard et al.ll2( 
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jol ISpringel et al.ll2005l; lAngulo et al.ll2008l: 


McBride et al.l 


2009; ICai et al. 20091: iTevssier et all 12009). 



Sanchez et all 120091 ; iPercival et all l2010l : iNorberg et al.l 



2011D . Such tools have been usually constructed following 



two methods, the first being ba sed on log-normal re aliza- 
tions of the matter density field (|Coles fc Joned [l99lD. and 
the second based on the use of N-body simulations. The log- 
normal catalogues can be produced in a shorter time period, 
and can be designed such that the final suite of mock ob- 
servations sh are the same cluste r ing properties a s that of 
real data (e.g lPercival et al.ll2004ICole et alj|2005l h On the 
other hand, N-body simulations have the advantage of prop- 
erly characterising the dynamical evolution and the growth 
of structures in the non-linear regime, providing a tool to 
test theoretical predictions for the abundance of dark matter 



In the context of galaxy clusters surveys, the combination 
of cosmological N-body simulations, a precise understand- 
ing of the survey selection function and a set of cluster-mass 
proxies is required for the construction of mock catalogues 
containing the relevant structural an dynamical properties of 
these objects, such as their masses, X-ray luminosities, spec- 
tral temperatures, fluxes and peculiar velocities, among oth- 
ers. This provides a control sample to analyse the possible 
systematic effects in the determination of scaling relations 
and/or cosmological parameters, together with the possibil- 
ity to assess the capability of future galaxy cluster surveys 
(e.g eROSITjQ) to push the constraints on the latter to 
higher precision (e.g lVikhlinin et al l (|2009h : iPillepich et al.l 

iSml)). 
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ing of galaxy clusters, prev ious analyses fe.g iRetziaff et all 
1 19981 ; ISchuecker et al.|[200lT ) used cosmological N-body sim- 
ulations in order to determine statistical properties such 
as covariance matrices associated to different measurements 
of abundance and clustering. On smaller scales, hydro-N- 
body simulations are widely used to explore the influence 
of non-thermal processes such as AGN a nd Supernova feed- 
back on the cluster scaling relations (e.g lEvrard et al.ll2002l ; 



Short et alj|2010h 



iBorgani et "al~l 120051 ; Istanek et al.ll201C 
The implementation of the baryonic physics required to 
model these astrophysical processes in large cosmological 
simulations is challenging and it would demand extremely 
large computational resources. This leads to the necessity to 
implement semi-analytic models, which in the case of galaxy 
clusters refers to the cluster scaling relations, such that the 
mock catalogues match the observed properties of these ob- 
jects. 

This is the second of a series of papers dedicated to the 
analysis of the new REFLEX II galaxy clust er catalogue. A 
first paper (jBalaguera-Antolfnez et all I^Oll! ). hereafter Pa- 
per I) was focused on the clustering analysis of this new sam- 
ple through the power spectrum. In that work, a set of mock 
catalogues was used to determine the covariance matrix of 
the cluster power spectrum and thereby assess the statisti- 
cal methods implemented in that analysis. In this paper we 
describe the construction of that suite of mock catalogues, 
which have been designed to reproduce basic properties of 
the REFLEX II sample (e.g survey geometry, selection func- 
tion, cluster abundance, etc.). A key ingredient in the con- 
struction of realistic cluster mock catalogues using N-body 
simulations is the assignment of luminosities to the dark 
matter haloes. 

Owing to the fact that the one-point statistics (e.g 
abundance) of galaxy clusters (as well as their clustering) 
are sensitive to the underlying mass M-X r ay luminosity Lx 
relation (hereafter M-Lx) (e.g lVoit]]2005l ), it is important 
to have a scaling relation consistent with both the observed 
cluster abundance and a halo mass-function in N-body sim- 
ulations. Accordingly, instead of assigning luminosities using 
scaling relations found in the literature (see for instance th e 
scaling relati o ns ob t ained bv | Reiprich fc Bohringer (1999); 
Stanek et al] (l2006h ; [Maughan I (|2007f ); iPratt etall lj2009li ; 
Mantz et al.l l|2010r )). we calibrate the M-Lx based on the 



X-ray luminosity function (XLF hereafter) of the REFLEX 
II sample. We discuss the main features of the scaling rela- 
tion, though it is not the scope of this paper to develop a 
deep analysis on this topic. 

We also briefly describe the two-point statistics of our 
mock catalogues (broadly described in Fourier space in Pa- 
per I) and show how these can be used to determine co- 
variance matrices in order to develop constraints on cos- 
mological parameters based on likelihood analysis (Sanchez 
et al., in preparation). The worth of our set of mock cata- 
logues is reflected in their physical and observational con- 
tent. These are translated to the covariance matrix of the 
cluster power spectrum (or correlation function) and its 
capability to properly account for non-linear evolution of 
the matter density field, dependencies of clustering strength 
with the underlying scaling relations and systematics ef- 
fects introduced by the survey selection function. In that 
sense, the set of REFLEX II mock catalogues introduced 



in this paper represents an improvement in t he analysis of 
this sample, compared to p revious works (e.g lCollins et all 
|2000|; ISchuecker et al.l|200lM based on the REFLEX sample 
l|Bohringer et al.ll200lf )~ 

The outline of this paper is as follows. In § [2] we briefly 
introduce the REFLEX II galaxy cluster catalogue and in 
§ 12.21 we describe the measurement and parameterisation 
of the REFLEX II XLF. In § [3] we introduce the N-body 
simulations and discuss the procedure we follow to calibrate 
the underlying M-L x relation (§|32]and §E3J. In §Hwe 
discuss some of the most relevant properties of the mock 
catalogues. Finally, we end with our main conclusions in 
§0 

We adopt a fiducial cosmological model based on a flat 
ACDM Universe with a matter energy density parameter 
of f2 ma t = 0.25, a dimensionless Hubble parameter h — 0.7 
(in units of 100kms _1 Mpc _1 ) and a rms of mass fluctua- 
tions of as = 0.773. Unless otherwise stated, redshifts, dis- 
tances, fluxes, luminosities, power spectra and mass func- 
tions are calculated with this set of cosmological parameters. 
Throughout this paper we always refer to the rest-frame X- 
ray luminosity in the ROSAT energy band 0.1 — 2.4 keV 
determined under the assumption of spherical overdensities 
(SO hereafter) characterised by A = 500 (i.e, the mean mat- 
ter density of the cluster is 500 times the critical density 
of the Universe). We will frequently refer to a set of re- 
alizations, namely, the N-body simulations of dark matter 
haloes. We warn the reader not to confuse these with the 
mock catalogues, which refers to the set of REFLEX Il-like 
catalogues. 



2 THE REFLEX II CATALOGUE 
2.1 The REFLEX II sample 

The REFLEX catalogue is bas ed on the ROSAT All-Sky 
Survey (RASS. ITruemper II1993T ). where galaxy cluster can- 
didates are detected as X-ray sources coinciding with galaxy 
concentrations in optical sky survey images (Bohringer et 
al., in preparation). The final identifications and redshift 
measurements come fro m compleme n tary follow-up obser- 
vations as described by iGuzzo et all (2009), yielding spec- 
troscopic redshifts for 860 clusters with flux limit of 1.8 x 
10 -12 ergs _1 cm -2 in the ROSAT energy band and approx- 
imately 6 per-cent incompleteness in redshift follow-up. The 
missing redshifts are currently being obtained though ob- 
servations at La Silla. While a detailed description of the 
construction of the catalogue and the selection function will 
be given in forthcoming papers (Bohringer et al. in prep.), 
Paper I presents a brief summary of the main steps followed 
in the derivation of the cluster parameters and the charac- 
terisation of the survey selection function. Here we briefly 
summarize the aspects that are more relevant to the con- 
struction of the mock catalogues. 

In our analysis we used the X-ray luminosities corrected 
for missing flux in the ROSAT energy band. The error in 
the measurement of the flux varies through the REFLEX II 
sample, depending on parameters such as the observed flux 
and the exposure time associated with each cluster. In order 
to keep the analysis simple, for the construction of the mock 
catalogues we adopted a normal distribution with a width 
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Figure 1. (a) The filled points with error bars correspond to 
the measurement of the REFLEX II X-ray luminosity function 
<E>(I/x)- The black dot-dashed line represents the best fitting 
extended-Schechter function, defined in Eq. JTJ, to this measure- 
ment. The red solid, blue-dotted and green dot-dashed histograms 
represent the mean luminosity distribution determined from the 
ensemble of 50 realizations of the L-BASICC II simulation using 
haloes with X-ray luminosities assigned according to i) a M-Lx 
relation given by Eq. Q (running index) with the set (a, b, c) cali- 
brated according to § 13,21 ii) a power-law M-Lx relation and iii) a 
power-law with mass-dependent scatter (§ 13,411 . respectively. The 
triangles represent the mean luminosity function from the final set 
of 100 mock catalogues. Panel (b) shows the difference between 
the measurements and the Schechter and extended Schechter fits 
in units of the standard deviation. 



cfl x = 0.2Lx (20 per cent error in the X-ray luminosity), 
which fairly characterizes the flux errors in the REFLEX II 
sample. 

The REFLEX II sensitivity map can be described by di- 
viding the surveyed area (13924 deg 2 in the southern hemi- 
sphere) into iVpix = 13952 pixels of approximately equal 
area of 1 deg 2 . Depending on its exposure time, the limiting 
luminosity L^ m (z) for each pixel was tabulated in the range 
0^2^ 0.8 assuming a nominal flux limit. This defines the 
REFLEX II selection function. 



2.2 The REFLEX II X-ray luminosity function 



As will be described in section [3721 a key ingredient for the 
construction of the mock catalogues is the X-ray luminosity 
function of the REFLEX II sample. We meas ured the XL F 
by means of the standard 1/Vmax estimator (|Felten|[l977l ). 
The resulting X-ray luminosity function, shown by the filled 
points with error bars in panel (a) of Fig. [TJ can be well 
described by what we call an extended-Schechter function: 



Here no determines the overall normalisation, a charac- 
terises the slope of the XLF in the low-luminosity regime 
and Li marks the transition from power-law behaviour at 
low luminosities to the exponential fall-off at high luminosi- 
ties. Th e function e„ (x) is the so-called g-exponential distri- 
bution (iTsallisI 120091 ) . defined as 



e q (x) 



(2) 



$(L)dL = no 



(1) 



The best-fitting parameters (no, a, L*, q) to the REFLEX II 
X-ray luminosity function were determined by implementing 
a M arkov-Chai n Monte Carlo (hereafter MCMC) algorithm 
(e.g |Neal|[l993T ). These are a = 1.54 ± 0.06, L* = (0.63 ± 
0.15) x 10 44 ergs/i -2 , n = (4.08T0.82) x 10~ 6 (Mpc /i" 1 )" 3 , 
and q = 1.31 ± 0.03, where the quoted errors denote the 
68 per cent confidence levels in each parameter. The best 
fitting extended-Schechter function corresponding to these 
parameter values is shown by the dot-dashed line in panel 
(a) of Fig. [T] Panel (b) of the same figure shows the differ- 
ence between the measured XLF and the best-fit averaged in 
the corresponding luminosity bins, in units of the standard 
deviation. For comparison, in the same panel we show the 
performance of the Schechter function^, showing that such 
parameterisation cannot fully characterize the REFLEX II 
XLF, especially at its bright tail. 

In order to carry out the statistical analysis using the 
information of the XLF, we assume a diagonal covariance 
matrix C* = o"|(Li)<5^ , with a variance given by the Pois- 
son standard deviation a^(Li) = $(Li) / \f~Ni, where Ni is 
the number of clusters in the i—th luminosity bin. This as- 
sumption will be also used within the analysis presented in 
§ 13.11 and represents a good description of the covariance 
matrix of the REFLEX II XLF, as will be shown in § H 



3 CONSTRUCTION OF THE REFLEX II 
MOCK CATALOGUES 

3.1 The N-body simulations 

Our mock catalogues are based on the Low resolution 
Baryon Acoustic Simulations at ICC (hereafter L-BASICC 
II) N-body simulations. A det ailed descript i on of these 
simulations c a n be found in lAngulo et all ((2008) and 
ISanchez et all (|2008l ). These correspond to 50 realizations 
of the same flat ACDM model characterised by a matter 
energy-density parameter S7 ma t = 0.237, a baryon energy- 
density f2ba = 0.046, a dimensionless Hubble parameter 
h = 0.73, a rms of mass fluctuations as = 0.773 and a 
scalar spectral index n s = 0.997. Each of the L-BASICC II 
simulations follows the dark matter distribution using 448 3 
particles over a comoving box of side 1. 34Gpc/i~ 1 , compa - 
rable to the Hubble Volume Simulation l|Evrard et alj|2002l ) . 
We used the halo samples at z = constructed by means of 
a friends-of-friends (hereafter FOF) algorithm with a linking 
lenght of b = 0.2. 

Panel (a) of Fig. [2] shows the mean halo mass function 



2 The best-fitting parameters of the Schechter function (q = 1) 
are a = 1.71 ± 0.07, L t = (1.51 ± 0.15) X 10 44 crgs/i~ 2 and 
n = (1.29 ±0.91) x 10- 6 (Mpch- 1 )- 3 
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n(M ) from the ensemble of simulations. Due to the resolu- 
tion of these simulations it is only possible to identify halos 
with masses M > 1.7 x 10 13 Mq/i -1 (corresponding to ten 
dark matter particles) , thus probing only the high- mass tail 
of the halo mass f unction, which i s very well described by the 
fitting formula of I Jenkins et al.l (|200ll ) (shown by the solid 
line in panel (a) of Fig.[2|. Note that the cosmological model 
of the L-BASICC II simulations is slightly different from our 
fiducial cosmology. We check that the results shown in this 
paper are not qualitatively modified due to these differences 
in cosmological parameters (see § 3.5). 

Even though the accuracy of the calibration of the pa- 
rameters in the different recipes for the mass function of dark 
mat ter haloes has been pushed to the 5 to 10 per cent leve l 
(e.g lJenkins et all 1200 ll ; I Warren et"al]|2006l ; ITinker II2008D . 
it has been shown that the inclusion of baryons in the N- 
body simulations in creases this uncertainty to ~ 15 per cent 
l|Stanek et al.|[2009l ). especially in the mass range probed by 
the L-BASICC II simulations. Such uncertainties must in 
principle be taken into account when comparing theoretical 
predictions with observations on cluster counts/abundances. 
Although we have not accounted for these uncertainties in 
the results presented in this paper, we checked that our con- 
clusions do not vary significantly when these are included by 
allowing up to a 20 per cent uncertainty in the normalisation 
of halo mass function. 

Similar ly, it has been shown (e.g ITinker I 120081 ; 

iLukic et all 120091 : iMore et al.l l201ll ) that the FOF halo 
finder systematically overestimates the masses of dark mat- 
ter haloes, compared to the masses of the same objects found 
within simulations at higher resolution. This can be a po- 
tential systematic effect that might affect the shape of the 
scaling relations linking the FOF halo masses to cluster ob- 
servables, in our case, the X- ray luminosity. Using the recent 
analysis of I More et al.l (|201 lh , we have verified that such cor- 
rections in the halo masses do not modify our results and 
conclusions (see § I3.4|l . 



3.2 X-ray luminosity assignment 

The first step in the construction of the REFLEX II mock 
catalogues from the N-body simulations is the assignment 
of X-ray luminosities to the dark matter halos. We do this 
following a M-Lx relation which we calibrate in order to re- 
produce the observed X-ray luminosity function of the RE- 
FLEX II sample. The information of the scaling relation M- 
Lx is encoded in the function p(Lx\M, z), specifying the 
probability of a cluster to have assigned a X-ray luminosity 
Lx conditional to its mass M at redshift 2. Given the abun- 
dance of dark matter haloes n(M,z), the X-ray luminosity 
function at a given redshift z is given by: 



$(L 



x,z = 



n(M, z)p(L x \M,z) dM. 



(3) 



In this way, given the halo abundance from the N-body sim- 
ulations, Eq. ([3| allows us to calibrate the M-Lx relation 
in order to reproduce the measured XLF. 

Under the assumption of virial and hydrostatic equi- 
librium, constant gas mass to total mass ratio and ther- 
mal bremmstralhung being the domin ant emissio n mecha- 
nism in the intra-cluster medium (see Voit 2005, for a re- 
view), the M-Lx relation should be well described by a 
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Figure 2. (a) Mean halo mass function n(M) (filled circles) of 
the L-BASICC II simulations at z = 0. Th e solid line shows the 
prediction based on the fit of I Jenkins et all l|200lh . The dotted 
and dashed lines show the effective mass function derived from a 
light-cone (following the REFLEX II sensitivity map (see § 13.40 ) 
for two different limiting X-ray luminosities, (b) M-Lx relation 
obtained from the abundance matching technique using the mass 
function at z = (dot-dashed line), and light-cone estimates using 
the FOF (solid line) and SO haloes with A = 500 (see § GQJ. 
(c) Slope (first derivative) of the M-Lx relation, (d) Curvature 
(second derivative) of the M-Lx relation. 



power-law, at least for the brightest galaxy c l usters (e.g 



Reiprich fc Bohringer! 1 19991: IStanek et aD 



2007 



Rvkoff et a l. 2008; Pratt et all I200E 



2006; Mauahan 



Vikhlinin et al 



20091 ). Furthermore, within the ROSAT energy band, the X- 
ray emissivity is weakly dependent on the temperature of the 
gas, specially for the hot clusters (T ^ 3 keV or Lx > L+). 
These facts can be shown to generate scaling relations of the 
form Lx oc M v ir and L^? 1 oc M*/ 3 for band and bolomet- 
ric luminosities respectively. Deviations from self-similar ex- 
pectations are associated with non-thermal processes such as 
AG N feedback from the growth of super-massive black-holes 



e.g iBorgani et al. Jl200a iKav et al. 1120071 : | Puchwein et al.l 



1200a IStanek et aLlboid lShort et aljboiol ). These feedback 
processes lead to a more significant gas depletion for low 
mass systems, reducing the X-ray luminosity and breaking 
the simple self-similarity as a function of mass. The evolu- 
tion o^he_j\4jj-ocjce^^ observations 
fe.g lVikhlinin et al.|[2009l : lReichert et alfeoill ) and found to 
also show deviations from the simple self-similar model. 

A simple non-parametric estimate of the underlying M- 
Lx relation can be obtained by impl ementing the so-calle d 
abundance matching technique (e.g iBehroozi et al.ll2010l ). 
in which we determine the X-ray luminosity Lx(M) as a 
function of the halo mass by matching the cumulative abun- 
dances <S>(L')dL' = n(M')dM'. This method is ap- 
plicable under the assumption that X-ray luminosity is a 
monotonically increasing function of the halo mass with no 
intrinsic scatter. The resulting M-Lx relation is shown by 
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the black dot-dashed line in panel (b) of Fig. [2] displaying 
deviations from a power-law. To obtai n this result we used 
the recipe for the halo mass function of lJenkins et all (|200lh 
computed for the cosmological model of the L-BASICC II 
simulations at z = 0. 

In order to quantify the inferred departure from a 
power-law M-Lx relation, the dot-dashed lines in panels 
(c) and (d) of Fig. [2] show the behaviour of its first and 
second derivatives, respectively, as a function of the halo 
mass. Within the luminosity range explored by the REFLEX 
II X-ray luminosity function, the slope of the M-Lx rela- 
tion displays a transition from a maximum value of approx- 
imately ~ 2 at ~ 5 x 10 13 Mg/i" 1 to ~ 1.4 at a mass scale of 
~ 10 16 Afgft _1 . Accordingly, the curvature displays positive 
values in the low halo mass regime M < 1Q 13 ' 5 MqK -1 . In the 
intermediate halo mass range 5 x 10 13 < M < 1O 15 M /i _1 , 
the curvature displays negative values and resume positive 
values for the high-mass haloes. Although these results are 
based on a simple halo abundance matching technique which 
ignores relevant physical effects, they suggest that a simple 
power-law relation M-Lx will not be enough to describe the 
observed XLF through Eq. . 

We calibrate the M-Lx relation that we use in the 
construction of our mock catalogues directly from the halo 
samples in the N-body simulations. We assume that X- 
ray luminosities are distributed around a mean scaling re- 
lation with a log-normal scatter in the luminosity (at a 
fixed mass) (J\ u l\m- In order to account for the deviation 
from a power-law as suggested by the abundance matching 
technique, we adopted a mass-dependent slope (or running- 
index scaling-law) to parameterise the mean M-Lx rela- 
tion in the whole luminosity range covered by the RE- 
FLEX II sample. Defining I = log 10 (Lx /10 44 ergs _1 h~ 2 ) 
and m = log 10 (M/10 14 Mo/i -1 ), our input mean scaling re- 
lation M-Lx has the form 



£(m) = a + bra + cm 



(4) 



where a determines the mean luminosity at 1O 14 M0/i -1 , 
b determines the slope at that mass-scale and the param- 
eter c characterises the curvature. We assigned X-ray lu- 
minosities to the dark matter haloes using Eq. Q, ap- 
plying a log-normal scatter around \nLx characterised by 
uin l\m = 0.26 (~ 30 per c ent scatter). This value was de- 
rived bv lStanek et al.l (|2006l ) based on cluster number counts 
and spatial clustering measured from a sub-sample of the 
REFLEX catalogue. We note however that such a value was 
determined under the assumption of a power-law sca li ng re - 
lation. Furthermore, in their analysis, Stanek et al.l l| 20061 ) 
assumed a fiducial cosmological model with fi ma t = 0.24 
and a power spectrum normalisation of erg = 0. 85, which 
is hig her than our fiducial value erg = 0.773). IReiprichl 
(2006 ) verified that the high intrinsic scatter measured by 
IStanek et ail <|200fJ) was related to the cosmological model 
they implemented, mainly, to the high value of erg. Using the 
cosmological param e ters derived fro m WMAP-3year data 
|Spergel et alj|2007f ). IReiprichl |2006f) found that a small in- 
trinsic scatter was still able to describe the observ ations 
(within the mass uncertainties). iMantz et al.l ((2010) used 
a fiducial cosmological model similar to ours, and found a 
larger scatter (~ 40 per cent), though it was substantially 
reduced (~ 6 per cent) when core-excised clusters were used 
in their analysis. We therefore emphasise that the param- 
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Figure 3. Distribution of the best-fitting parameters (a, b, c) de- 
termined from the 50 realizations of the L-BASICC II simulations 
within the full luminosity range as well as in the faint and bright 
tail of the XLF. Solid line distributions show the results assum- 
ing a M-Lx relation given by Eq. Ht . Filled histograms show the 
distribution of parameters associated with a power-law (c = 0). 



eters of the M-Lx relation obtained with a fixed scatter 
and from a FOF halo mass function are only meant to re- 
produce the cluster abundance within the N-body simulations 
and z = 0. 

On top of the intrinsic scatter, we introduced the er- 
rors in the X-ray luminosity through a second log-normal 
scatter with o = 0.2. These two scatters can be added in 
quadrature, such that the total M-Lx relation preserves a 
log-normal scatter given by a = 0.328 (39 per cent scatter). 
After luminosities were assigned to each halo we determined 
the luminosity functions within the 50 simulation boxes. We 
implemented a MCMC algorithm to determine the set of 
parameters (a, b, c) in Eq. Q by comparing the obtained lu- 
minosity distributions with the measured REFLEX II XLF, 
assuming a Poisson-noise standard deviation, as pointed out 
in §[221 

As a consistency check, we have verified that replacing 
the measured XLF by the best fit extended-Schechter func- 
tion of Eq. fTJ yields a set of parameters (a, b, c) different by 
less than 5 per cent. 



3.3 The inferred M-Lx relation 

Following the procedure described in § 13.21 we obtained best- 
fitting values of the parameters (a, b, c) from each of the 50 
L-BASICC II realizations. The distributions of these values 
are shown in the upper panels of Fig. [3] and present narrow 
peaks around the values a — — 1.36±0.03, b — 1.88±0.05 and 
c = —0.29 ± 0.04 (where the errors represent the variance 
from the ensemble of simulations). These values, together 
with the effective scatter <r = 0.328, were used in the final 
construction of the mock catalogues described in § 13.51 and 
used in Paper I. The constraints on c indicate a ~ 7a devi- 
ation from the power-law M-Lx relation (i.e. c = 0). The 
solid line in Fig. [4] shows the mean M-Lx relation obtained 
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Figure 4. M-Lx relation of the REFLEX II mock catalogues. 
The solid line represents the mean M-Lx relation calibrated from 
<J>(Lx). Filled circles correspond to individual halos in one of 
the L-BASICC II simulations with X-ray luminosities assigned 
as described in § 13. 5^ . Open circles correspond to mock clusters 
in one member of the ensemble of 100 REFLEX II mock cata- 
logue. The dashed line represents the M-Lx relation determined 
by the abundance matching technique (§ 13,51) . The horizontal line 
marks the value L*. The inner plot shows the distribution of 
the luminosities with respect to their mean value. The solid line 
represents a Gaussian distribution with dispersion a centred at 
ln(L/L) = (d In V/d In L)a 2 ~ 0.14, where V(L) is the volume as 
a function of the X-ray luminosity (sec § I4.1|l . 



with these set of parameters. The points correspond to in- 
dividual halos in one of the L-BASICC II realizations with 
X-ray luminosities assigned following our M-Lx relation, 
while open circles represent clusters in one of our final RE- 
FLEX II mock catalogues (i.e, once the REFLEX II selection 
function is applied as described in § 13. 5[) . As shown by the 
solid histogram in Fig. [1] the luminosity distribution gener- 
ated through this calibration provides an excellent match to 
the measured X-ray luminosity function of the REFLEX II 
sample. 

The M-Lx relation obtained through this calibration 
procedure is in agreement with that derived from the simple 
abundance matching technique described in § 12.21 (shown by 
the dashed line in Fig.Q}, especially for M < 5 x lO 14 M /i _1 
which corresponds to a luminosity Lx ~ L+. At higher 
masses, the difference between the two is clear. Due to the 
exponential fall-off of the halo mass function at such high 
mass scales, the obtained X-ray luminosity distribution of 
the halos is very sensitive to the assumed scatter in the M- 
Lx relation. 

As a consistency check, we have verified that the values 
(a, b, c) obtained from the resulting distribution within the 
L-BASICC II simulations are in good agreement wit h those 
obtained when the fitting formulae of I Jenkins et al.l (|200ll ) 
is used to compare the measured &(Lx) against the pre- 
diction from Eq. ©. The results obtained in this case are 
a = -1.29 ± 0.04, b = 1.95 ± 0.09 and c = -0.31 ± 0.06, 
where the errors denote the la confidence intervals. 



A power-law scaling relation could still be sufficient to 
reproduce the abundance of high luminosity (Lx > L*) 
galaxy clusters. In order to verify this, we repeated the pro- 
cedure described in § 13.21 to calibrated the M-Lx relation 
in the faint (Lx ^ £*) and bright (Lx > L*) tails of the 
luminosity function. The results are shown in the second 
and third rows of Fig(3] The set of parameters calibrated 
in the first interval (Lx < £*) are marginally compatible 
with a power-law (for instance ~ 6 per cent of the realiza- 
tions displayed curvatures in the range |cj ^ 0.05), whereas 
the analysis derived from the bright tail is more consistent 
with such scaling relation (~ 25 percent of the realizations 
had |c| ^ 0.05). We have checked that by pushing the limit- 
ing luminosity up to higher values, the resulting distribution 
of the curvature is well peaked at the value c ~ 0, though 
the distribution becomes wider due to the low number of 
objects with such luminosities. Fig. [3] also shows the distri- 
bution of the parameters associated with a power-law M-Lx 
relation, using the full range as well as the faint and bright 
end tails of the luminosity distribution. The evident shifts 
in the peaking values of the slope b in these two luminos- 
ity ranges shows the need for a mass-dependent slope, or in 
general, for deviations from a power-law M-Lx relation in 
case the full luminosity range is used in the analysis. 



3.4 On the shape of the M-Lx scaling relation 

As discussed in § [2] the determination of the X-ray lumi- 
nosities of the REFLEX II sample assumed, among other 
things, a fiducial set of cosmological parameters and a def- 
inition of halo masses, namely, SO masses with A = 500. 
In this subsection we briefly explore the possible impact of 
the differences between these fiducial values and the corre- 
sponding quantities characterising the L-BASICC II N-body 
simulations, on the recovered values of the parameters of the 
M-Lx relation. 



• Resolution effects. I t has been shown ( e.g lWarren et all 
120061 ; iLukic et all 120091 ; iMore et al.l 1201 lh that the finite 



number of particles defining a dark matter halo is overes- 
timated by the FOF finder algorithm, when compared to 
the same masses obtained with simulations of higher resolu- 
tion. I n particular, the empirical correction of lWarren et al.l 
(2006), based on the number of particles defining a dark 
matter halo, was designed to match the abundance of low 
resolution to that of higher resolution simulations. This cor- 
rection leads to a 40 per cent mass overestimation at the 
resolution limit of the L-BASICC II simulations, reducing 
to a few percent at 1O 15 M0 hT 1 . Based on perc olation the- 
ory, th is systematic effect has been modeled bv lMore et al.l 
l|201ll ). who showed that the FOF masses are overestimated 
by ~ 30 per cent at the resolution of the L-BASICC II simu- 
lations, and by 15 per cent at a mass scale of 10 15 Mq/i -1 (S. 
More, private co mmunication) , repres enting a higher correc- 
tion than that of lWarren et al.l |2006t ). In order to check the 
impact of this systematic effect on our results, we corrected 
the halo masses within the L-BASICC II simulations and 
re-calibrated the scaling relation given by Eq. (4[ on a num- 
ber of realizations. The results show a 15 per cent increment 
in the normalization a, with its corresponding decrease the 
slope b (around 10 per cent) and a shift of the curvature c 
towards positive values (around 15 per cent), with respect to 
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FOF L-BASICC II Simulation -1.36 ± 0.03 1.88 ± 0.05 -0.29 ± 0.04 

FOF L-BASICC II Simulation: masses corrected -1.10 ± 0.02 1.70 ± 0.04 -0.23 ± 0.04 

FOF fjenkins et al. 20011 -1.29 ±0.04 1.95 ± 0.09 -0.31 ± 0.06 

SO A = 500 (Tinker 2008) full L x range -0.68 ± 0.06 1.65 ± 0.03 -0.38 ± 0.04 

SO A = 500 (Tinker 20081 L > 1.6L* -0.69 ± 0.07 1.45 ± 0.26 -0.15 ±0.22 

SO A = 500 (Tinker 20081 L > 1.6L*, power-law -0.64 ± 0.03 1.27 ± 0.05 



Table 1. Parameters of the M-Lx relation in Eq. Q obtained by calibration against the observed REF LEX II XL F: i) using the 
L-BASICC II simu lations, ii) using the L-BASICC II simulations with FOF masse s correc t ed as iMore et al. | fcOllll , iii) using the 
I Jenkins et al.| | |200lf ) fitting formula for the halo mass function, and iv) & v) using thelTinkerl feOOSf) mass function for SO masses with 
A = 500 in two luminosity intervals. In order to compare with the analysis of |]VIant2^r*alT j2010T) , the last row vi) shows the results 
obtained by fitting a power-law in the range L > 1.6L*. 



our original estimations. In particular, the shift in the curva- 
ture corresponds to approximately la, excluding the value 
c = with a statistical significance of ~ 4<r. We can con- 
clude that the correction for the finite number of particles 
in some extent diminishes the departure from a power-law 
scaling relation, although not completely. 

• Cosmological parameters. We have tested the impact 
of the cosmological parameters in the calibration of the 
scaling relation, by constraining the set (a, 6, c) using the 
IJenkins et all (|200ll ) mass function at 2 = computed for 
different values of f2 ma t in the range 0.15 — 0.35, keeping 
both as fixed to our fiducial value or letting it vary as 
a 8 = 0.45f2~° t 3 (e.g lEvrard et al.ll2002T l. In both cases, while 
the amplitude a (the slope 6) decreases (increases) with re- 
spect to fimat, the curvature c remains approximately con- 
stant, with a 2a significance. 

• Redshift evolution. The redshift evolution of the mean 
matter density of the Universe c an be imprinted in the am - 
plitude of the M-L x relation (e.g lBorgani ll200ll ; IVoitll2005t ). 
The parameterisation given by Eq. Q can then be thought 
of as accounting for such evolution through the conversion 
A(z) — > 10 a M cm , in case the true underlying M-L x relation 
is given by a power-law Lx = A(z)M b . On the other hand, 
even without a redshift dependent amplitude (i.e, A(z) con- 
stant), the function p(Lx\M, z) could still evolve with red - 
shift through its intrinsic scatter (e.g lFabiian et all 1201 lT ) . 
For the purposes of building mock catalogues from a snap- 
shot of the simulations, such dependency could be translated 
to a mass-dependent scatter. We explored this possibility by 
calibrating a power-law M-Lx relation with a scatter given 
by a{M) = a lnL | M (Af/10 14 M Q /!- 1 ) 7 . The resulting X-ray 
luminosity distribution generated by this parameterisation 
(with best-fitting values a ~ —1.2, b ~ 1.8 and 7 ~ —2.3 
with ~ 50 per cent error in 7) is shown by the dot-dashed 
histogram in Fig. [TJ Such parameterisation represents a mild 
improvement with respect to the power-law or running in- 
dex M-Lx relation in the bright end (L > L*) of the XLF, 
while working as good as the power-law for L < L+ clusters. 

• Light-cone effects. In order to take into account the red- 
shift evolution of the halo abundance, the XLF measured 
within a past light-cone must be compared with the vol- 
ume average of Eq. Under the assumption of a redshift- 
independent M-Lx relation, we can write the prediction for 
the luminosity functions in the light-cone as 



where n e g (M, 2 ma x) is the halo mass function averaged over 
the volume of the light-cone defined by the sample: 



n eS (M, z„ 



V(z n 



maxj Jy(z max ) 



dU 

n(M,z)—dz, 
dz 



(6) 



§{L X )= n cff (M,z max )p(Lx|M)dM, 



(•») 



with z m ax = Zmax(ix) as the maximum redshift allowed 
for a cluster with luminosity Lx in order to be detected, 
and U(z m a X ) is the comoving volume at that redshift. This 
information is encoded in the REFLEX II sensitivity map 
described in § 12.11 As an example, we show the behaviour 
of n e ff (M, Zmax) in panel (a) of Fig. [2l based on the halo 
mass function for FOF and SO masses. In order to apply 
the match abundance technique to derive an estimate of 
the shape of the M-Lx relation using Eq. ((6|, we carried 
out an iterative process in which an initial power-law scal- 
ing relation was assumed, £ = —1.0 + 1.5 m, from which 
n c ff (M, z max (Af )) could be determined in order to match 
the abundances, leading to a new estimation of the M-Lx 
relation, which is used as an input for the next iteration. The 
resulting M-Lx relations (convergence is achieved after few 
iterations) derived from the FOF and SO mass functions are 
shown in panel (b) of Fig. [2] with its first and second deriva- 
tives in panels (c) and (d) respectively. The recovered scaling 
relation still displays departures from a power-law, although 
the curvature is slightly smaller than that derived from the 
simulation snapshots (dot-dashed line in Fig. [2}. The de- 
viation from a power-law is more evident for FOF masses 
than for SO masses. This might indicate that the deviations 
from a power-law scaling relation observed by matching the 
abundances at redshift z = are not significantly accounted 
for by the evolution of the halo-abundance within the light 
cone denned by the REFLEX II sample, at least in the case 
of vanishing intrinsic scatter. 

• Definition of halo mass. We calibrated the M- Lx rela- 
tion g iven by Eq. @ using the fitting formulae of iTinkerl 
(2008) for SO haloes with A = 500. The mass function was 
evaluated at the median redshift of the sample z « 0.09. 
Fixing the intrinsic scatter to our fiducial value we obtained 
a = -0.68 ± 0.06, b = 1.65 ± 0.03 and c = -0.38 ± 0.04. 
The value for the curvature excludes the power-law scal- 
ing relation with ~ 7a confidence. Similarly, for clusters 
with Lx > Li, (with a median redshift of z w 0.14) 
we found a = -0.73 ± 0.02, b = 1.68 ± 0.16, and c = 
—0.34 ± 0.14, excluding the power-law behavior to ~ 2a 
significance. We calibrated a power-law scaling relation in 
the range Lx > 1.6L* ~ 10 44 erg s _1 h~ 2 , similar to that 
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Figure 5. Slice of lOOMpc/i -1 width through the centre of the 
halo distribution in one of the L-BASICC II simulations. The 
solid lines indicate the geometric configuration used to extract 
two independent mock catalogues with z max = 0.22 from each 
realization. 
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Figure 6. Properties of the mock catalogues: (a) right-ascension, 
(b) declination, (c) redshift (d) luminosity (e) luminosity relative 
to the true underlying M-Lx relation and (f) mass distribution. 
The shaded areas show the standard deviation in the distribu- 
tion of each property derived from the set of mock catalogues 
and compared to the corresponding property in the REFLEX II 
sample (solid line). In panel (f) the solid line shows a Gaussian 
distribution centred at 8a 2 (see § 14.11 1 



used by iMantz et all (|201Ch , obtaining a = -0.64 ± 0.05 
and b = 1.27 ± 0.05, in good agreement with that analysis. 
In the same range, the fit to Eq. Q yields a curvature of 
c = —0.15 ± 0.23, compatible with a power-law. 

The tests presented in this section show that the deviation 
from a power-law scaling relation in the full X-ray luminos- 
ity range probed by the REFLEX II XLF is likely not due to 
differences in the mass definition, the cosmological param- 
eters or evolution along the light cone. This is nevertheless 
not conclusive, in the sense that a full analysis including the 
evolution of the scaling relation, its intrinsic scatter and the 
uncertainties in the halo abundance is still required in or- 
der to draw stronger conclusions about the cluster scaling 
relations. 

Table [TJ summarizes the values of the parameters (o, b, c) 
obtained from the various analysis presented in this section. 

3.5 The construction of the mock catalogues 

We now describe the procedure followed to construct the 
REFLEX II mock catalogues. First, we note that in Pa- 
per I it was shown that the effec tive volume for clustering 
analyses fe.g lFeldman et aI1ll994h probed by the REFLEX 
II sample does not increase substantially beyond redshift 
2 « 0.22. As the primary goal of the set of mock catalogues 
is to determine covariance matrices for the statistical anal- 
ysis of the large-scale structure probed by the REFLEX II 
sample, we set this value as the maximum redshift of the 
mock catalogues. This in turn allowed us to construct 100 
independent mocks catalogues out of the 50 L-BASICC II re- 
alizations. Mock catalogues with higher maximum redshifts 
can in principle be constructed, at the cost of having less 
independent realizations. 

The geometric arrangement of the mock catalogues is 



depicted in Fig. [5] where a slice cut of 100 ft _1 Mpc through 
the center of one of the L-BASICC II realizations is shown. 
Given such geometric arrangement, we assigned a cosmo- 
logical redshift z cos to each halo by inverting the redshift- 
comoving distance relation r(z cos ) = J Zcos dz'H~ 1 (z') for 
the cosmology of the simulations. The final redshifts of the 
halos were obtained taking into account the distortions in- 
troduced by their peculiar velocities v as z = (1 + z C os)(l + 
Zpcc) — 1, where z pcc = v/Hq. For the z — output of the 
L-BASICC II simulations, the distribution of peculiar veloc- 
ities displays a mean value of v ~ 350kms _1 with a veloc- 
ity dispersion of a « 180kms _1 , which in terms of redshift 
co rresponds to Az ~ 6 x 10~ 4 . As can be seen in table 2 
of Guzzo et al l l|2009h . the redshift distortions induced by 
spectroscopic redshift uncertainties in the REFLEX II sam- 
ple are significantly smaller than those induced by peculiar 
velocities. Therefore we do not introduce such effect in the 
estimate of the redshift of our mock clusters. 

We assigned X-ray luminosities to the dark matter 
haloes according to the M-Lx relation described in § 13.31 
These halos are then observed through the REFLEX II sen- 
sitivity map, rejecting those haloes with luminosities below 
the limiting luminosity L^ m at a given location. The open 
circles in Fig. [4] correspond to individual objects in one of 
our final REFLEX II mock catalogues which, with a min- 
imum mass M ~ 3.5 x 10 13 /i _1 Mq, corresponding to a 20 
particle dark matter halo, covers a mass range well within 
the resolution limit of the simulations. 



4 PROPERTIES OF THE REFLEX II MOCK 
CATALOGUES 

In this section we discuss some of the most relevant proper- 
ties of the REFLEX II mock catalogues. 
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Figure 7. Covariance matrix of the XLF. Panel (a) shows the 
comparison between the standard deviation of the XLF derived 
from the set of mock catalogues, and the Poisson-noise prediction 
determined from the number of clusters in each bin of X-ray lu- 
minosity. Panel (b) shows the correlation coefficients r%j of the 
XLF. 



4.1 Basic properties 

The underlying M-Lx relation applied to construct the RE- 
FLEX II mock catalogues (and used in the clustering anal- 
ysis of Paper I) is characterised by Eq. (j4]), with the set of 
parameters shown in § 13.31 The performance of the REFLEX 
II sky-mask is shown in panels (a) and (b) of Fig.[S]by means 
of the distribution in equatorial coordinates of the mock cat- 
alogues, compared to that of the REFLEX II sample, shown 
by the solid line. The Milky-way band is well described in 
panel (a) by the two minimum distribution values at a ~ 100 
and a ~ 250 deg (see also fig. 1 of Paper I). Panels (c) and 
(d) show respectively the redshift and X-ray luminosity dis- 
tributions, displaying good agreement with that of the real 
REFLEX II sample. The mass distribution within the mock 
catalogues is shown in panel (f) of Fig. [B] This implies an 
average halo mass of M ~ 4 x 10 14 /i _1 Af s for the REFLEX 
II clusters. 

The application of the REFLEX II selection function to 
the N-body simulations, introduces in our set of mock cata- 
logues a selection effect inherent to all flux-limited samples, 
namely, the fact that objects with increasing luminosities 
probe larger cosmological volumes. In combination with the 
intrinsic scatter in the M-Lx relation this introduces the so- 
called Malmquist bias which, if unaccounted for, can lead to 
biased estimates of the underlying scaling relation, especially 
when a threshold in luminosity is imposed to the sample. 
Our set of mock catalogues can then be used as a control 
sample to test different methods to obtain scaling relations 
taking into account the selection effects mentioned above. 
At a given mass M the Malmquist bias causes the mean 



observed M-Lx relation (L(M)) to deviate fro m the under- 
lying M-Lx relation L(M ) according to (e. g IStanek et al.l 
12009 ; IVikhlinin et~ai1l2009l ; iPratt et al.ll2009r i 

(lnLx) =\nL x (M) + 8a 2 , (7) 

where it is assumed that the volume probed by a cluster 
of luminosity Lx scales as V(Lx) oc L s x (5 = for a vol- 
ume limited sample) together with a log-normal distribution 
function p(Lx\M, z) with scatter a. Under the assumption 
of a non-evolving M-Lx relation, the Malmquist bias can 
be translated into an overestimation of the slope and an un- 
derestimation of the ampli tude in the M-Lx relation (e.g 
iReiprich fc Bohringer|[l999h . In panel (e) of Fig. © we show 
the distribution of luminosities around the mean M-Lx re- 
lation given by Eq. @, for the 100 mock catalogues. The 
dashed line represents a log-normal distribution centred at 
\n(L/L) = 8a 2 « 0.14, where we have used 8 = 1.38, which 
gives a good description of the maximum volume as a func- 
tion of the X-ray luminosity. In terms of luminosities, this 
bias is translated into a factor ~ 1.15. 

We finally study the behavior of the covariance matrix 
of the REFLEX II XLF. In panel (a) of Fig. [3 we compare 
the standard deviation obtained from the final set of mock 
catalogues with that derived from a Poisson noise. The good 
agreement between these two estimates implies that the in- 
homogeneous spatial distribution of X-ray galaxy clusters 
(i.e, the sampling variance) has a minor imp act on the de - 
termination of the variance of the XLF (e.g ISmith II2012T ) 
given the comoving volume and luminosity range probed by 
the REFLEX II sample. Furthermore, we verified that the 
off-diagonal elements of the covariance matrix of the XLF 
are negligible compared to the variance. In this regard, panel 
(b) of Fig. [7] shows the correlation coefficients defined by 

\rtj\ = \Cfj |/ >J C*Cjj- determined from the suite of mock 
catalogues, showing in general weak correlations between 
different luminosity bins, ranging from nj ~ 10~ 3 at low 
luminosities to m ~ 10~ 10 at high luminosities A more de- 
tailed analysis of the REFLEX II XLF will be presented in 
a forthcoming paper (Bohringer et al., in preparation). 

4.2 Clustering properties 

The power spectra of the REFLEX II mock catalogues were 
analysed in Paper I. These measurements were obtained 
following the method describe d in Paper I, b ased on the 
minimal- variance estimator of iFeldman et all (|l994 ). This 
study showed that the mean power spectra from our ensem- 
ble of mocks catalogues for different limiting luminosities is 
in excellent agreement with those of the REFLEX II sam- 
ple on intermediate and large scales (0.02 < fc/(/iMpc -1 ) < 
0.3). While a detailed analysis of the clustering in configura- 
tion space will be presented in a forthcoming paper (Sanchez 
et al., in preparation), also in this case we find good agree- 
ment between the clustering signal of the REFLEX II sam- 
ple and that of the mock catalogues. This allows us to use 
our mock catalogues for its main application in the context 
of large scale structure, that is, to obtain estimates of the 
covariance matrices of these measurements. 

Paper I presented an analysis of the bin-averaged covari- 
ance matrix C(ki,kj) of the cluster power spectrum. This 
covariance matrix showed highly correlated modes (with 
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Figure 8. Variance of the power spectrum (panel a) and the 
correlation function (panel b) measured from the mock catalogues 
(solid lines) compare d with the prediction from the estimator of 
iFeldman et alj l ll994t) (dashed lines), for two values of minimum 
X-ray luminosity. 



correlation coefficients r%j > 0.9) on intermediate (k ~ 
0.08ft Mpc -1 ) and small scales (k ~ 0.3ft Mpc -1 ), mostly 
due to mode coupling ind uced by the survey window func- 
tion. I^ldma^^^al] ( 19941 ) derived an approximated expres- 
sion for the variance of the spherically averaged power spec- 
trum under the assumption of Gaussian-distributed power 
around the true underlying value, given by: 



(8) 



a 2 P (k) _ 2 

P(fc) 2 ~ VfcVeff(fc)' 

where Vk ~ 47rfc 2 <5fc/(2-7r) 3 is the volume of a spherical shell 
of width 8k and V e ff(fc) is the effecti ve (coherence) volum e 
probed by the survey at a scale k (e.g lFeldman et al]|l994l ). 
In its simpler form (for a volume limited sample), this con- 
tains the contribution from cosmic variance and shot-noise. 
The prediction of Eq. (JSj) can easily be transla ted into the 
variance of the two-point correlation fu nction as (|Cohnl2006l : 
ISmith et al.ll2008l ; ISanchez et al.ll2008l ): 

4(r) = ^J o dkk 2 j (kr) 2 a 2 P (k), (9) 

where jo(y) is the zero-th order spherical Bessel function. 

Panel (a) of Fig. [8] compares the variance of the clus- 
ter power spectrum derived from the set of mock catalogues 
using sub-samples with minimum luminosities L™ ln = 0.049 
and L™ m = 0.588 x 10 44 ergs _1 ft" 2 against the prediction 
of Eq. |(S} . Similarly, panel (b) of the same figure shows the 
behaviour of the variance of the cluster correlation function 
for the same values of minimum X-ray luminosity. Despite 
its simplicity, Eq. (|5J| gives a good account of the variances 
inferred from our mock catalogues. However, significant dif- 
ferences exist, especially for the small-scale a^, which high- 
light the importance of the ensemble of mock catalogues as 
a tool to obtain accurate estimations of covariance matrices 
which take into account the full effect of non-linearities. 

Our ensemble of mock catalogues allows us to test an- 
other common assumption when comparing clustering mea- 
surements against the predictions of cosmological models, 
that of a Gaussian likelihood function. By studying the dis- 
tribution function of the values of P(k) and £(s) obtained 
from the mocks it is possible to reconstruct the shape of 
the likelihood function and find deviations from the sim- 
ple Gaussian assumption. As the covariance matrices for 
these statistics are not diagonal, this analysis is simplified by 
working on their "de-correlated" versions Xi and Q defined 
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Figure 9. Probability distribution function (PDF) of the values 
of the power spectra (panel a) and correlation functions (panel 
b) in our ensemble of mock catalogues for two values of minimum 
X-ray luminosity. The dashed lines represent a Gaussian PDF. 



\ i = (Mp l ) ij P j , 



(10) 

(ii) 



where Mp and are the transformation matrices of the 
basis where the covariance matrices Cp and are diag- 
onal (given by their eigenvectors) . Fig. [9] shows the prob- 
ability distribution function of AA = (Aj — Ai)/<T» and 
A(" = (fj — Ci)/fi f° r two limiting luminosities for the power 
spectrum (panel a) and correlation function (panel b) . These 
are well described by the prediction for a Gaussian distri- 
bution, shown by the dashed line, indicating that deviations 
from the Gaussian case are very small. 

Fig. [TJJ shows the distribution of the chi-squared val- 
ues x 2 associated with the cluster power spectrum and the 
correlation function within the 100 mock catalogues, com- 
puted with respect to their mean values within the ensemble. 
The ob served distributions are well d escribed by a \ 2 distri- 
bution l|Abramowitz fc Stegun Ill972l ) shown by the dashed 
lines, computed using the number of modes in each scale 
interval shown in the figure. This implies that the distribu- 
tion of these quantities at a fixed scale is likely to follow 
a Gaussian distribution. We have t ested this by perf orm- 
ing a Kolmogorov-Smirnov test (e.g lPress et al.ll2002l ). as- 
suming as a null hypothesis a Gaussian distribution for the 
quantity (P(fej) — P(ki))/P(ki), with zero mean and dis- 
persion a(ki)/P(ki). The test shows that on scales 0.02 < 
ki/(h Mpc" 1 ) < 0.2 the distribution is compatible with the 
null-hypothesis to a > 88 per cent level. We can therefore 
conclude that the resulting covariance matrix of the RE- 
FLEX II power spectrum can be used in an error analysis 
based on the assumption of Gaussian likelihoods. Similar 
conclusions can be obtained with respect to the clustering 
in configuration space. 
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Figure 10. Distribution of the x 2 values associated to the power 
spectrum (panel (a)) and the correlation function (panel (b)) of 
the 100 mock catalogues computed with respect to the mean of 
the ensemble. The dashed line represents the x 2 distribution for 
the corresponding number of modes in each interval of scales or 
wavenumbers. 



5 CONCLUSIONS 

In this paper we describe the procedure to construct a suite 
of galaxy cluster mock catalogues based on the X-ray lumi- 
nosity function $(Lx) of the REFLEX II sample and dark 
matter halo samples identified in the L-BASICC II N-body 
simulations. 

The construction of our mock catalogues relies on the 
scaling relation between halo masses and X-ray luminosity in 
the ROSAT energy band. Given the fact that the abundance 
of galaxy clusters is highly sensitive to the underlying scaling 
relation, we calibrate our mass proxy based on the X-ray 
luminosity function of the REFLEX II sample and the halo- 
abundance of the N-body simulations. We find that in the 
luminosity range probed by the REFLEX II sample, the 
observed $(Lx) demands a M-Lx relation that deviates 
from a power-law. 

We qualitatively discarded such behaviour as a by- 
product of the differences between the features of the N- 
body simulations and the assumptions under which the RE- 
FLEX II catalogue was built, namely, the definition of halo 
masses and the cosmological parameters. Similarly, the devi- 
ations from a power-law are likely not associated with possi- 
ble light-cone effects (at least under the assumption of negli- 
gible intrinsic scatter). However, we find that a power-law is 
still sufficient to generate the observed abundance for bright 
clusters ( Lx ^ 10 44 erg s h~ 2 ), in agreement w ith recent re- 
sults (e.g lPratt et al.ll2009l ; liviantz et alj|20ich . This might 
suggest that the observed break in the slope of the M-Lx 
relation is a direct consequence of the physical properties of 
the galaxy clusters within the REFLEX II sample. A thor- 
ough analysis taking into account both the evolution of the 
halo mass function and the evolution of the full M-Lx re- 
lation is required in order to draw more solid conclusions 



about the shape of the scaling relation as demanded from 
the observed abundance. 

Our set of mock catalogues reproduce the observed 
abundance of the REFLEX II sample as well as the ob- 
served two-point clustering statistics of this data-set. There- 
fore, they provide an excellent tool to the test the statistical 
methods applied to the real data and to perform error anal- 
ysis, especially through the covariance matrix of the cluster 
power spectra or the correlation function. The covariance 
matrix derived from this suite of catalogues properly ac- 
counts for physical effects such as the non-linear evolution 
of clustering, its connection to the underlying scaling rela- 
tions, as well as the effect of the survey window function. 
This represents a major advantage of our set of mock cata- 
logues against theoretical estimations of covariance matrices 
or log-normal catalogues. Our results also show that the dis- 
tribution function of the values of P(k) and £(s) obtained 
from the individual mocks are well described by the pre- 
diction for a Gaussian distribution. This indicates that it 
is possible to assume a Gaussian likelihood function when 
comparing measurements of these statistics with the predic- 
tions from cosmological models. 

In the coming years, the next generation of galaxy clus- 
ter surveys will probe much larger volumes than any present- 
day sample. For example, the eROSITA cluster survey is 
expected to contain approxi mately 10 5 objects (w ith ~ 10 3 
with measured redshifts) (e.g lPillepich et al]|201ll ) spanning 
a wide range of masses over a volume of ~ 10 2 Gpc 3 . The 
construction of the mock catalogues necessary for the anal- 
ysis of the large-scale clustering of this sample will demand 
a great effort. This task will require a suite of extremely 
large N-body simulations with a unique combination of vol- 
ume and resolution to construct mock catalogues in the form 
of light-cones up to a redshift ~ 1.5. The implementation of 
baryonic physics in such simulations would be even more am- 
bitious and these mock catalogues will need to rely on more 
indirect methods like the one discussed here. Although the 
volume of the REFLEX II sample allowed us to implement a 
few simplifications such as the use of a single snapshot of the 
simulations, the scheme we have implemented to calibrate 
the M-Lx relation can be extended to light-cones spanning 
wider redshift ranges, providing a simple technique to con- 
struct mock catalogues out of pure dark matter simulations 
even for these large samples. 
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